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Based on the Wilemski-Fixman approach (J. Chem. Phys. 60, 866 (1974)) we 
showed that for a flexible chain in 9 solvent hydrodynamic interaction treated with 
an pre- aver aging approximation makes ring closing faster if the chain is not very short. 
Only for a very short chain the ring closing is slower with hydrodynamic interaction 
on. We have also shown that the ring closing time for a chain with hydrodynamic 
interaction in 6 solvent scales with the chain length (N) as iV L527 , in good agreement 
with previous renormalization group calculation based prediction by Freidman et al. 
(Phys. Rev. A. 40, 5950 (1989)). 

I. INTRODUCTION 



Dynamics of loop formation in long chain molecules has been a subject of immense interest 
to experimentalists jl-4| and theoreticians js 17 1 for over a decade. Recently the dynamics of 
loop formation has found greater importance because of its relevance in biophysics. Advances 
in single molecule techniques made it possible to monitor the kinetics of loop formation 
involving biomolecules at the single molecule level 18|, LL9|. Loop formation is a prime step 



in protein [20] and RNA folding 



21]. Also a measure of the intrinsic flexibility of DNA can 



18]. 



be found from the rate at which it undergoes cyclization 

Loop formation in polymers is essentially a many body problem as a polymer is made 
of many connected segments and hence an exact analytical solution is impossible. All the 
theories of loop closing dynamics are approximate |s , |s|| . Although many theoretical and 
simulation attempts have been made to investigate the effect of flexibility [17j, solvent quality 



13| on the loop closing dynamics, not much theoretical investigation has been performed to 



shed light on the effects of hydrodynamic interaction on loop closing dynamics other than the 
renormalization group calculation by Freidman et al [22]. Hydrodynamic interaction which 



2 



is essentially nonlocal in space has been shown to have profound effects on the dynamics of 



long chain molecules 



23 



24j . The rate of translocation of polymer through a nano-pore has 



24 



25). Recently 



been shown to be greatly affected in presence of hydrodynamic interaction 

it has been theoretically shown [23] that the breakage rate of stretched polymer tethered to 

soft bond get enhanced in presence of hydrodynamic interaction. 

In this paper we investigate the effect of hydrodynamic interaction on the ring closing 

dynamics of a flexible chain in 9 solvent. It is well known that the ring closing time (r) 
or a flexible chain without excluded volume and hydrodynamic interaction (Rouse chain 
26]) scales with the chain length (N) as N 2 in the Wilemski-Fixman (WF) [5] theoretical 



framework. Here we use Zimm model [26j] for the polymer which actually takes care of 
the hydrodynamic interaction at the simplest possible level. Our calculation shows that 
hydrodynamic interaction profoundly affect the rate of loop formation and also has a different 
scaling relation with the length of the polymer N. The rest of the paper is arranged as 
follows. In Sec. II the Wilmeski-Fixman (WF) theory for the chain closure jljj is briefly 
discussed. WF formalism gives a prescription to calculate the ring closing time as an integral 
over a sink-sink time correlation function. Sec. Ill deals with a brief description of the radial 
delta function sink which is used later in the calculation of the ring closing time. Time 
correlation formalism for the flexible chain with and without hydrodynamic interaction is 
discussed in Sec. IV. Sec. V presents the results and VI is devoted to conclusions. 

II. THEORY OF CHAIN CLOSURE 

In an stochastic environment the dynamics of a single polymer chain having reactive 
end-groups is modeled by the following Smoluchowski equation in WF theory 

gP({ * },t) = LP({R}, t) - A;5({R})P({R}, t) (1) 

Here P({R},t) is the distribution function for the chain that it has the conformation 

{R} = i?i,i?2 ; Rn at time t where Ri denotes the position of the ith monomer in the 

chain of n monomers. S(R) is called the sink function which actually models the reaction 
between the ends and thus usually is a function of end to end vector. L is a differential 
operator, defined as 



3 



d 



d dU 
+ 



P({R},t) 



(2) 



dR, dR 

Here Do is the diffusion coefficient of the chain defines as the inverse of the friction 
coefficient per unit length and U is the potential energy of the chain. Wilemski and Fixman 
then derived an approximate expression for the mean first passage time from Eq. (CQ). This 
mean first passage time is actually the loop closing time for the chain. The expression for 
this loop closing time reads 



T 



dt 



C(t) 



1 



v C(oo) 

Where C(t) is the sink-sink correlation function defined as 



(3) 



C(t) 



dR / dR S(R)G(R,t|R ,0)S(R )P(Ro 



(4) 



In the above expression G(R, t|R ,0) is the Greens function or the conditional probability 
that a chain with end-to-end distance Rq at time t = has the end-to-end distance R at 
time t; P(Rq) is the equilibrium distribution of the end to end distance since the chain was 



in equilibrium at time t — 0. S(R) is the sink function |27H29| which depends only on the 
separation between the chain ends. We would like to comment that Eq. ([3]) is not exact and 
only valid in the limit of infinite sink strength, k — > oo as this closing time is nothing but 



the mean first passage time 



30 



3l|. 



Now it is obvious that the knowledge of the Greens function, G(R, t\R Q , 0) is prerequisite 
to calculate the sink-sink correlation function C(t) and hence the closing time r. In case of 
a flexible chain the Greens function and the end-to-end probability distribution functions 
are known and Gaussian. 

For a flexible chain the Greens function is given by 



3/2 



G(R,t|R ,0) 
Where 



x exp 



3(R-0(t)R o ) 2 
2( J R%(i-0 2 W) 



(5) 



(R(t).R(O)) 



eq 



(iP> 



(6) 



eq 
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is the normalized end-to-end vector correlation function for the chain. The above ensem- 
ble average is taken over the initial equilibrium distribution for end-to-end vector P(R ). 



Similarly end-to-end equilibrium distribution for the flexible chain at time t = is given 



by 



15 



P(R ) 



3/2 



exp 



2L 2 



(7) 



2ttL 2 J 

With the above Gaussian functions the sink-sink correlation function can be written as 
a radial double integral. 



C(t) 



2nL 2 



(i-m 3/2 J ° 



4nR 2 S(R)dR / 4nR 2 S(R )dR x 



exp 



3 (R 2 + R 2 



sinh [m(t)RRo)/(L 2 (l - 2 (t)))] 



2L 2 (1 _ j (34>(t)RRo)/(L 2 (l - 2 (t))) 



(8) 



The above integral can be evaluated analytically for some specific choice of the sink 
functions. With a radial delta function sink the above integral can be evaluated analytically. 

III. THE RADIAL DELTA FUNCTION SINK 

For our case study we choose a radial delta function sink [9], S(R) = 5(R — a). Since in 
this case the integration over R and Rq can be carried out analytically, the looping time can 
be expressed in a closed form as follows. 



T 



» dt [ exp [-(2x o 2 (t))/(l - 2 (t))] sinh [(2x <f>(t)) / (1 - 2 (t))] 

U \ AU\ /l JL2/ 



(9) 



2x o 0(t)Jl-0 2 (t) 



with 



x 



3^ 
2L 2 



Obviously if the end-to-end vector correlation function <p{t) is known the closing time can 
be calculated by carrying out the integration over time. Here we calculate the looping time 
given by the above expression for a flexible polymer in presence and absence of hydrodynamic 
interaction and make a comparative study. Throughout the paper it is assumed that the 



5 



chain was in equilibrium at t = and only at t = + the hydrodynamic interaction is turned 
on (for the Zimm chain). Thus the equilibrium end-to-end distribution for the Rouse as well 
as for the Zimm chain is given by (Eq.flJJ)). The last integration over time is not analytical 
so has to be carried out numerically. 



IV. FLEXIBLE POLYMER WITHOUT AND WITH HYDRODYNAMIC 

INTERACTIONS 

The simplest dynamical description of a flexible chain or polymer in solution is given 
by Rouse Model 



26 



32j. This model does not take into account of excluded volume and 
hydrodynamic interaction and has been the basis of dynamics of dilute polymer solution 
and has been successfully used in the context of ring closing in polymers [lOj, breathing 
dynamics in DNA \%\ . polymer translocation [34] etc. In the continuum limit the governing 
equation of motion for the position vector R(n, t) is given by 

Here n is a continuous variable, k = 3k ^ T where b is the length of a monomer in the discrete 
representation of the Rouse polymer and f(n,t) is the random force with the moments 



(f(n,t)) = 

(f a (n, h)f p {m, h)) = 2(k B T5(n - m)5 a pS(t 1 - t 2 ) (11) 

It is straightforward to show that the Rouse normal modes (X p ) obey the following 
equation 

C-Qf = - k vX P + f P {t) (12) 
Where C = N(, ( p = 2N( for p = 0, 1, 2, ... and k p = 2ir 2 kp 2 /N = ^|fV for p = 0, 1, 2, .... 
Here f p (t) is the random force satisfying (f pa ) = and {f pa (t)f q p(ti)) = 2S pq S a pC p kBT5(t — 
fi). 

Then for a Rouse polymer it can be shown that the normalized end-to-end vector time 
correlation function defined in (jSJ) is given by 
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= < R y. = V 4- 2 e-^ = V -^e-*^" (13) 

Where n = 4 I Stt and c^ ouse = v -. 

The simplest possible model for the dynamics of a flexible chain with the hydro dynamic 
interaction but without the excluded volume interaction is known as the Zimm model 26 



In 9 condition the continuum limit equation of motion for the position vector R(n, t) in the 
Zimm model is given by 

Here H mn is the mobility matrix and is generally a nonlinear function of R(n, t) — R(m, t) 
and the above equation is quite difficult to handle. To simplify this analysis Zimm introduced 
a pre-averaging approximation, which replaces H mn by its equilibrium average, (H mn ) . It 
can be shown that with this pre-averaging approximation 26] becomes a linear function of 
R(n,t). 

BH(n,t) . . / d 2 R(m,t) „. A 

= X>(»- ™) [k ^ ) + /K «) J (15) 

where h(n — m) decreases slowly as h{n — m) oc \{n — m)\^ 1 ^ 2 . Thus in Zimm model 
interaction among the segments is not localized. This is how a Zimm chain in different from 
a Rouse chain. 

It is possible to show 26] that for the Zimm chain in 9 solvent, the normal modes 
(Xp imm ' d ) can be shown to obey an equation which has the same structure as the Rouse 
normal modes 



C P P m = -k p X z p ^ e + f p (t) (16) 

Where C p = C\f^f L for P = Q > 1 , 2 > ■■■ and K = 27r 2 kp 2 /N = ^|fV for p = 0, 1, 2, .... Here 
f p (t) is the random force satisfying (f pa ) = and (f pa (t)f qf3 (ti)) = 25 pq 5 al3 ( p k B T5(t -ti). 

In case of 9 solvent condition the normalized end-to end vector time correlation function 
(Eq.(jSJ)) for the Zimm chain is given by 
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Figure 1: <p(t) against time (t). The values of other parameters used are: N = 1,6 = 1,£ = 
l,fceT = 1. In the inset the short time (t) behavior of r is shown. Note up to a short time 
range <p(t) for the chain with hydrodynamic interaction (Zimm chain in 8 solvent) decays slowly as 
compared to a chain without hydrodynamic interaction (Rouse chain). 



= (R(*)-R(0)>e« = y 8 _^ /Tf w = ^ 8 e _ twfmm , 

Wed 



_Zimm,6» _ C^ 372 ^ a „A ,.,Zimm,8 _ P 3/2 

Zimm, 6 



Where T-^ rmm > w — ' flnf ] ,,,2? 

Where r x - 6v ^3/2 fcsT and u p 

In principle it is straightforward to calculate the closing time with a delta function sink 
(Eq.Q) as it involves evaluating <f>(t) and carrying out an integration over time. We evaluate 
4>(t) for both the cases exactly by carrying out the sums denned in (Eq. (TT3]) ) and (Eq.flTj 



V. RESULTS 



The closing time (r) with the radial delta function sink is calculated for the Rouse chain 
as well as for the Zimm chain in 9 solvent. Calculation of r involves evaluating 0(t) for the 
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Figure 2: (f)(t) against time (t). The values of other parameters used are: N = 20, b = 1,( = 
l,k B T=l. 

respective chains and then putting it in Eq. ([9]) followed by an integration over time which 
has to be carried out numerically. The end-to-end vector correlation function defined in Eq. 

is actually a sum over the decays of all the odd modes and at a given time t the bulk of the 
contribution to <p(t) comes from the lower normal modes. Each normal mode decays with an 
effective rate constant, Up OUSe = p 2 jr\ for the Rouse chain (Eq.(TT3l)) and in case of a Zimm 
chain in 6 solvent it is u;f mm ' e = p 3 / 2 / T f imm > e (Eq.jnj)). Now had it been n = T f mm ' 8 all 
the odd modes of the Rouse chain hence the end-to-end vector correlation function would 
have decayed faster resulting a faster ring closure compared to that for a Zimm chain in 
6 solvent. In reality T\ ^ T ^ tmm ^ e anc l the ratio { J j^ j tmm ^ /u)p OUSe has a strong length (N) 

/ i Zimm, 9 \ , — , — 

dependence, u Zvmm,e j u Rouse = >- = M^^N^. Only when, ^/Wpp « ^ the 

normal modes of the Zimm chain in 9 solvent decay faster. For a very short chain or in case 
of higher normal modes the condition y/N/p « can be achieved and then only the 
normal modes for the Zimm chain (in 6 solvent) will decay slowly compared to the normal 
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Figure 3: r against the chain length(iV). The values of other parameters used are: b = 1, £ = 
1, feT = 1, a = 0.01. In the inset r for a short range of the chain length, iV shown. Note for a very 
short chain the r for the chain with hydro dynamic interaction (Zimm chain in 9 solvent) is higher 
meaning a slower ring closure. 

modes of the Rouse chain. Thus for a short Zimm chain (in 9 solvent) initial decay of <p(t) is 
slow compared to a Rouse chain. This is shown in the inset of Fig. 1. On the other hand for 
a long Zimm chain all the normal modes decay faster than the Rouse chain normal modes. 
This naturally results a slowly decaying <j)(t) for a Rouse chain as is shown in Fig. 2. 

In Fig. 3 log-log plot of the closing time (r) is plotted against the chain length N for 
both the chains. The slope for the Rouse chain has been found out to be 2.003 very close 
to previously found values based on WF theory calculation 
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For the Zimm chain 



(in 9 solvent) the slope is 1.527. Thus hydrodynamic interaction makes ring closure faster. 
Also notice for a very short chain the ring closure is faster for a Rouse chain as shown in 
the inset of Fig. 3. This actually suggests that the proportionality factor or the frequency 
factor between the ring closing rate, r and N a , where a = 1.527 with the hydrodynamic 
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interaction on and a = 2.003 without the hydrodynamic interaction, too has a length (N) 
dependence. 

VI. CONCLUSIONS 

In this paper we show that the hydrodynamic interaction makes ring closure faster within 
WF [5] framework if the chain is not very short. For a very short chain ring closure can 
be slower in presence of hydrodynamic interaction. The ring closing time (r) for the Zimm 
chain (in 9 solvent) scales as iV 1,527 . This is actually in very good agreement with the 
renormalization group calculation prediction of Friedman et al. jz^]. They showed that 
for a Zimm chain not too short the closing time scales as iV 3 / 2 . But surprisingly N 3 / 2 
behavior can also be seen even without hydrodynamic interaction. For example, harmonic 
chain approximation within WF framework leads to similar scaling Jf|. Also in the model of 
Szabo, Schulten and Schulten (SSS) [8fl, in which the difficult problem of the polymeric chain 
having many degrees of freedom is replaced with a single particle diffusing in a potential of 



mean force predicts TV 3 / 2 scaling 



br a chain without hydrodynamic interaction. Portman 



14] pointed out that SSS theory 



8| actually gives a lower bound to the chain closing time 
and WF theory [5|] gives an upper bound. So a faster ring closure in SSS formalism as 
compared to WF theory prediction is expected. Experimentally the closure rate k for a 
10 — 20 residue polypeptides of the alanine-gylcine-glutamine trimer have been found out 
to be vary as N~ 3 ^ 2 for large N [3]. Although in accordance with SSS theory [8], this faster 
ring closing as compared to a flexible chain without hydrodynamic interaction might also 
be due to hydrodynamic interaction as found in our theoretical analysis here. 

Thus just by looking at the length dependence of the ring closing time it is rather impos- 
sible to comment on the microscopic basis of ring closure. Similar scalings can arise due to 
completely different reasons. As a future problem it will be interesting to see what happens 
to the closing time for a chain with hydrodynamic interaction as different sink functions are 
used and also in the SSS formalism |8J. It is expected since the SSS formalism gives a lower 
bound to the ring closing rate [14j, ring closing with the hydrodynamic interaction work in 
the SSS framework will make ring closing even faster as compared to in the WF framework. 
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